arXiv:1508.01000vl [math.OC] 5 Aug 2015 


Noname manuscript No. 

(will be inserted by the editor) 


Uniform Quadratic Optimization and Extensions 

Shu Wang • Yong Xia 


Received: date / Accepted: date 


Abstract The uniform quadratic optimizatin problem (UQ) is a nonconvex 
quadratic constrained quadratic programming (QCQP) sharing the same Hes¬ 
sian matrix. Based on the second-order cone programming (SOCP) relaxation, 
we establish a new sufficient condition to guarantee strong duality for (UQ) 
and then extend it to (QCQP), which not only covers several well-known re¬ 
sults in literature but also partially gives answers to a few open questions. For 
convex constrained nonconvex (UQ), we propose an improved approximation 
algorithm based on (SOCP). Our approximation bound is dimensional inde¬ 
pendent. As an application, we establish the first approximation bound for the 
problem of finding the Chebyshev center of the intersection of several balls. 
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1 Introduction 

We study the nonconvex uniform quadratic optimization problem: 

(UQ) max ier fo(x) 

s. t. k < fi(x ) <Ui, i = 1,... ,p, 
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where — oo < U < Ui < +oo and the quadratic functions fi(x) are defined by 

fi(x) = x t Qx + 2 bfx + di , i = 0, 1 ,... ,p, 

with Q being a real symmetric matrix, 6 * € R" and di € R for i = 0, 1 ,... ,p. 
Throughout this paper, we denote by (UQ+) the widely used special case of 
(UQ) where Q >- 0 (a notation standing for that Q is positive definite). 

The problem (UQ+) was first introduced by Beck [3i to find the smallest 
ball enclosing a given intersection of balls, or equivalently, find the Chebysliev 
center of the intersection of given balls, which can be reformulated as the 
following nrin-max problem: 


(CC) minnrax||x 

2 x€f2 


z\\ 2 = min 

Z 


{ max 


2 x T z) + || 2 || 2 


( 1 ) 


where Q = {x € R" : ||a; — a,;|| 2 < rf, * = 1,... ,p}, || • || denotes the Euclidean 
norm, and the inner maximization problem is a special case of (UQ) with Q 
being the identity matrix. Moreover, we note that any NP-liard binary integer 
linear programming problem 


T 

max c x 

s. t. af x <bi, i = 1,..., to, 

* e {o, i} n , 

where c € R n , a* € R" and 6,gl for i = 1,,m, can be reformulated as the 
following special case of (UQ+): 

max x x + (c — e) x 

s. t. x T x + (a* — e) T x <bi , i = 1 ,..., to , 

0 < x T x — e T x < 0 , 

0 < x T x + (ej — e) T x <1, j = 1,..., n, 

where e € R” is the vector of all ones and ej £ R™ is the j-th column of the 
identity matrix. Since the 0-1 Knapsack problem is already NP-hard, we see 
that the problem (UQ+) remains NP-hard even when 

p > n + 2. (2) 

Lagrangian duality plays a fundamental role in optimization, especially in 
the quadratic constrained quadratic programming (QCQP) problem. It pro¬ 
vides an upper bound of the primal maximization problem, i.e., weak duality 
holds ■ It is a common sense that strong duality holds for convex optimiza¬ 
tion under Slater condition, that is, the difference between the optimal values 
of the primal and dual problems is zero j7]. Though, it fails for the general 
non-convex programs, strong duality may hold for some nonconvex (QCQP), 
see for example, pilFHT(ni^ll^Hll^llMHi?H] and references therein. 
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For the problem (UQ+), Beck [3] showed that the strong duality holds as 
long as 

p < n — 1. (3) 

Replacing the inner maximization problem of © with its Lagrangian dual 
minimization problem, Beck [3] established the convex relaxation of ©. Under 
the above Assumption (|5jl. the convex relaxation is tight and hence the min- 
max problem © is globally solved. Later, Beck [?] further relaxed the above 
assumption m to 

p<n, (4) 

and then proved that the strong duality holds for the problem (UQ+) if for 
each i = 1,... ,p, exactly one of the following three cases occurs: 


(i) k = -oo, (ii) k = Ui, (iii) it* = +oo. 


(5) 


Consequently, the min-max problem © is globally solved by Beck’s convex 
relaxation under the weaker assumption ©. However, when p > n, Beck’s 
convex relaxation is no longer tight, and, to our knowledge, the quality of 
Beck’s bound is still unknown. 

In this paper, we will make a thorough study on the problem (UQ). Our 
main contributions can be divided into the following three parts. 

(i) We pay a revisit to the problem (UQ+) without assuming exactly one 
of the three cases © holds (see Section 2). Differently from the semidefinite 
programming (SDP) reformulations we establish the following second- 

order cone programming (SOCP) relaxation: 


(SOCP) 


max x£ Kr. t + 2 bgX + do 

s. t. U < t + 2 bf x + di < Ui, i = 1,... ,p, 


Q 2 x 

t-i 

2 


< 


t + 1 


( 6 ) 

(7) 


and show that it provides an upper bound as tight as the primal SDP relax¬ 
ation. Using a simple and easy-to-understand proof, we show the equivalence 
between (UQ+) and (SOCP) under the following assumption 


either rank [b \,..., b p ] < n — 1 or p = n, (8) 

which is slightly weaker than Assumption ©. Besides, under the additional 
primal Slater assumption, we establish the strong duality for (UQ+) and its 
Lagrangian dual. 

(ii)We extend the above SOCP reformulation approach from (UQ+) to 
the general non-convex (QCQP) (see Section 3). More precisely, we establish 
a new sufficient condition under which (QCQP) is hidden convex, i.e., it is 
equivalent to a convex programming problem. As applications, we show that 
our new sufficient condition not only generalizes a few existing results for 
special (QCQP) but also partially gives answers to a few open questions in 
literature. 








4 


Shu Wang, Yong Xia 


Our first corollary is that the trust region subproblem (TRS) jT2j enjoys 
the strong duality, which is a well-known result mmm- Moreover, our new 
convex reformulation brings a new look at the hidden convexity of (TRS). Our 
sufficient condition for the strong duality of the extended trust region problem 
with linear inequality constraints coincides with Hsia and Sheu’s condition m 
and improves Jeyakumar and Li’s condition m- As a further extension, we 
establish a new sufficient condition for the hidden convexity of variants of the 
extended trust region subproblem with not only linear constraints but also 
two-sided ball constraints [9]. 

For the weighted maximin dispersion problem mm, our sufficient condi¬ 
tion answers the open question [TS] under what condition the corresponding 
semidefinite programming relaxation is tight. 

As another application, we obtain the first sufficient condition to guarantee 
the hidden convexity of the general problem (UQ) without assuming Q >~ 0. 

Finally, we consider the extended p-regularized subproblem (p-RS) [TH251 
mm with additional linear constraints. A class of polynomial solvable cases 
of the extended (p-RS) has been studied in [TB] where p = 4 and the number 
of linear constraints is fixed as a constant. The case p / 4 remains unknown 
|18| . Our sufficient condition identifies a class of polynomially solved cases of 
the extended (p-RS) with any p > 2. 

(iii)We propose an improved approximation algorithm for convex-constrained 
(UQ+), i.e., U — — oo and iq < +00 for i = 1,... ,p (see Section 4). Actually, 
for the nonconvex quadratic optimization problem with ellipsoid constraints: 

(ECQP) max Ig fn f(x) = x T Ax + 2 b T x 

s. t. || F k x + g k \\ 2 < 1 , k = 1 ,... ,p, 

where A £ M. nxn is symmetric, b £ R™, F k £ R'’ , ' xn , g k £ iF* and ||p fc || < 1 
for k = 1,... ,p, Tseng m proposed an approximation algorithm based on 
the SDP relaxation to find a feasible solution x of (ECQP) in polynomial time 
such that 

/(*) > ■ ^ SDP )> ( 0 ) 

where 7 := maxfc = i i ... iP ||p fc || and v(-) denotes the optimal value of problem (•). 
Very recently, Hsia et al. [20] improved the approximation bound ([9]) to 

f(x)>(^^-) -t;(SDP), (10) 

Vvr + 7/ 

where r = min 

holds for the convex-constrained (UQ+), as it is a special case of (ECQP). 

In this paper, based on our second-order cone programming relaxation, we 
propose a new approximation algorithm, which finds a feasible solution x in 
polynomial time such that 

/(l)i (7f^) 2 -'' (SOCP) ' 


V8p+17-3 

2 


, n + 1 >. Trivially, this approximation ratio m 


(ii) 
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where 7 := maxi = i p . ^ We notice that the approximation 

V^- d *+IIQ^MI 2 

bound cut) greatly improves (HOD as we can show that v (SOCP) — v (SDP). 
It should be noted that, though this new approximation bound relies on the 
input data, it is independent of the numbers p and n. 

(iv)We extend the above new-developed approximation analysis for the 
convex-constrained (UQ+) to the problem of finding the smallest ball enclos¬ 
ing a given intersection of balls, i.e., the min-max problem (jT]) (see Section 
5). Replacing the inner maximization problem with its Lagrangian dual re¬ 
laxation, Beck proposed an efficient convex quadratic relaxation [§], which 
globally solves m when p < n. But it could fail to find the global minimizer 
for the hard case p > n. To our knowledge, the quality of Beck’s convex 
quadratic relaxation remain unknown as well as the quality of the feasible so¬ 
lution returned by Beck’s approach. In this paper, we succeed in establishing 
the first approximation analysis. We remark that our approximation bound 
dependents only on the distribution of the given balls rather than p and n. 

Throughout the paper, v(-) stands for the optimal value of problem (•). Let 
K" denote the n-dimensional vector space. Let I r the r x r identity matrix. 
For a matrix Q , denote by Af{Q) and 7 Z(Q) the null and range space of Q, 
respectively. A m i n (Q) and A max (Q) stand for the smallest and largest eigen¬ 
values of Q, respectively. Q >- (^)O means that Q is positive (semi)definite. 
For Q y 0, denote by Q 2 the square root of Q, i.e., Q = Q^Q?. Moreover, 
if Q >~ 0, is the inverse of Q 2 . The inner product of two matrices A 
and B is denoted by A • B = Tt(AB t ) = Y^i =1 Sj=i II ’ II denotes the 

Euclidean norm, i.e., ||x|| = Vx T x. For a set 17 C KL ra , int(l7) denotes the set 
of all the interior points in 17. The linear subspace generated by {b±,, b p } 
is denoted by span{&i,..., b p }. dim{M} denotes the dimension of the linear 
subspace M. 


2 An SOCP Relaxation and Strong Duality for (UQ_|_) 

In this section, we present a second-order cone programming relaxation (SOCP) 
for (UQ+). It is a common sense that solving an SOCP is much easier than 
SDP. Theoretically, though it is equivalent to the SDP relaxation, (SOCP) 
not only provides a simple and easy way to characterize the strong duality of 
(UQ+), but also implies a slightly weaker sufficient condition for the strong 
duality. 


2.1 SDP Relaxations 

We first consider the Lagrangian dual problem of (UQ+). Introducing p free 
Lagrange multipliers Ai,..., A p yields the Lagrangian function of (UQ+): 
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L(x, A) = f 0 (x) + ^2 A+ (ui - fi(x)) + ^2 A, ; {fi(x) - k) 

i= 1 i= 1 

( V \ / P \ T P 

1 - ^ Ai J x t Qx + 2 I 6 0 - ^2 x i b i ) % - ''22 x i d i 

P 

+ Y Ui — \ + do, 


i=1 


i=l 


i=l 


where A^ = max{Ai,0}, A, ; = — min{A,;, 0}. and note that Aj = Xf — X i for 
i = 1,... ,p. Then the Lagrangian dual problem of (UQ_|_) can be written as 


(D) inf i d(A) := sup L(x, A) 

AeRr ( xeR n 

By using Shor’s relaxation scheme m, the dual problem (D) can be reformu¬ 
lated as the dual semidefmite programming (SDP) relaxation for (UQ+): 


(D-SDP) inf r + ^(A^ Ui - k) - ^ A id* + d 0 

i= 1 i=l 

Q bo 
bo d 0 

Xi e R, i = l,...,p. 


3. t. J2 A * 


2=1 

r G 


Q bi 
bT di 


0 0 
0 T 


h 0, 


The primal SDP relaxation for (UQ+) can be directly obtained by lifting 
x g K™ to Y = xx T £ K. raxn and then relaxing Y = xx T to Y Y xx T (see [3]): 


(P-SDPi) 


max 


f Q bo 
Uo do 




s. t. U < 


f Q bi 
\bf di 


Y x\ 
x T 1 ) 




<Ui, i = ,,p, 


( 12 ) 


It is similar to verify that (D-SDP) is also the conic dual problem of (P- 
SDPi). Since <3^0, there is a strictly feasible solution for (D-SDP). If we 
further assume that the Slater condition holds for (UQ+), then (P-SDPi) also 
has a strictly feasible solution. Therefore, strong duality holds for (D-SDP) 
and (P-SDPi), i.e., v(D-SDP)= u(P-SDPi) (see, e.g., [5Tj). 

By introducing the invertible transformation y = Q 1 x and Ci = Q~^bi for 
* = 0,1,... ,p, we reformulate (UQ+) as 

(UQ 0 ) maxjgp y T y + 2 c^y + d 0 

s. t. U-di< y T y + 2 cjy <Ui~d z , i = 1,... ,p. 
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Lifting y £ R" to (^b !/) £ jjdn+i)x(n+i) anc [ then relaxing y T y to 

z £ R, we obtain a new SDP relaxation for (UQo) [3]: 


(P-SDP 2 


max 


Onxn Co 

Cn 1 


In y 

T 

y z 


+ do 


s. t. li — di < 


Onxn C, 

cT i 


In y 

T 

y z 


<Ui-di, i = 1 ,... ,p, 


In y 

T 

y z 


y o. 


Beck jniHI use d (P-SDPi) and (P-SDP 2 ) to prove the strong duality under 
Assumptions © and ©, respectively. 


2.2 An SOCP Relaxation 


To simplify the presentation, we first assume di = 0 for i = 0,1,... ,p. Other¬ 
wise, for % = 1 ,... ,p, we replace Ui and U with ut — di and h — di, respectively. 
We also assume Q = I n as it has been made in (UQo). Then, (UQ+) is reduced 
to the following formulation: 

(U) max x£ Rr, f 0 (x) = x T x + b^x 

s. t. U < x T x + 2 bf x < Ui, i = 1 ,... ,p. 


The next result is the main theorem of the section. 


Theorem 1 Under the Assumption m, the problem (U) is equivalent to the 
following SOCP relaxation: 


(S) 


max x6 r t + b^x 

s. t. k <t + 2bfx < Ui, i = 1 ,... ,p, 



(13) 

(14) 


in the sense that x* globally solves (U) if and only if (x*,t*) := (x*,x* T x*) 
globally solves (S). 


Proof. We first notice that the constraint (fT41) is equivalent to 

x T x < t. 


(15) 


Suppose v(S) = Too, it follows from (fT5l) that for any unbounded feasible 
solution of (S), say (x,t), we have t —>■ +oo. Since 

t + 2bjx >t — 2 \JbJbiVx T x = t — 2||& i ||v / t —> +oo, 

we have Ui = +oo, for i = 1,... ,p. Thus, we obtain u(U) = +oo. Now we can 
assume u(S) < +00 and let (x*,t*) be an optimal solution of (S). If ( x*) T x* = 
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t*, then x* is a feasible solution of (U) and v(S) = fo(x*) < u(U). On the other 
hand, as (S) is a relaxation of (U), we always have u(S) > u(U). Therefore, 
it holds that u(S) = u(U), which completes the proof. In the remaining, we 
assume ( x*) T x * < t* and consider the following two cases: 

(i) Suppose rank [&i,..., b p ] < n — 1. Let m be the largest number such that 
there are m indices ii ,..., i m £ { 1 ,... ,p} satisfying that rank [b ^,..., 6 j m ] = 
m and the ik -th constraint in m is active at (x*,t*) for k = 1 
According to the dehnitions, we have m < n — 1 and 

t* + 2bJ k x* =S ik , k = 1 ,..., to, 

where Si k = U k or Ui k for k = 1,... ,m. Then, there is a nonzero vector 
b G M" such that bf k b = 0 for k = 1,..., to. Since 

x(e) := x* + eb 

remains feasible for (S) when |e| is sufficiently small, we have b T bo = 0. It 
follows that x(e) remains optimal for small |e|. Define 

E = {e G R | k <t* + 2 bjx{e) < m, i = 1,... ,p}, 

J = {j 6 { 1 ,... ,p} \ {* 1 ,... ,imj | rank [b ^,..., 6 j m , bj\ =m+ 1 , bJb^O}. 

If E is unbounded, there is an e such that x(e) t x{£) = t* and hence 
completes the proof. Otherwise, it is sufficient to assume J / 0 and E is 
bounded, since J = 0 implies the unboundedness of E. In this case, there is 
a jo G J and £o such that x(eo) is feasible for (S) and the jo-th constraint 
in m becomes active at (x*,t*), we obtain a contradiction against the 
definition of to. 

(ii) Suppose p = n. According to Case (i), it is sufficient to consider the case 
rank [b \, •.., b p ] = n and all the n constraints (1131) are active at (x*,t*), 
i.e., 

Bx* = S — t*e, (16) 

where S = (Si) E R n with Si = U or m, e G is the vectors of all ones, 
the coefficient matrix B = [26i, 262 ,..., 2 b n ] T is nonsingular. Define 

x(t) = B~ 1 (5 — te). 


Then (S) is equivalently reduced to the following one-dimensional opti¬ 
mization problem: 

(S') max t + b£x(t) = max {b^B~ 1 5 + t(l - 6 ^B“ 1 e)} , 

t^xit) 1 x(t) g {t)<.0 

where g(t) = (S — te) T B~ T B _1 (S — te) — t is a strictly convex quadratic 
function. It follows from the facts 

x(t*) = x *, g(t*) = x* T x* — t* < 0, lim g(t) = +00 

£—>• + 00 
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that the equation g(t) = 0 has two roots t\ < t 2 and 
{t | g(t) < 0} = {t | ti < t < t 2 \. 

Since the objective function of (S') is linear, either t\ or t 2 is an optimal 
solution of (S'). But both satisfy t = x(t) T x(t). The proof is complete. 

□ 

The following example illustrates that Assumption © seems to be tight 
and can not be expected for an improvement. 

Example 1 Consider the following instance of (U) 

max x 2 

s. t. 1 < x 2 + 2x < 3, 

— 1 < x 2 — 2x < 3, 

where n = 1, p = 2, and rank[bi,b 2 ] = 1- Assumption © fails to hold. It is 
easy to verify that the optimal value is 1. However, the SOCP relaxation (S) 


max t 

s. t. 1 < t + 2x < 3, 


-l<t-2x<3, 



has an optimal value 3. 

The strong duality for (UQ+) was first established in [3] under Assumption 
© and then strengthened in [T] by assuming that © holds and exactly one of 
the three cases © occurs. At the end of this section, we show that, following 
Theorem [2] the second-order cone programming relaxation © is as tight as 
the primal SDP relaxation (fl2l) and the strong duality immediately holds for 
(UQ+) under the slightly weaker Assumption ©. 

Theorem 2 Suppose the Slater condition holds for (UQ , ) and u(UQ,) < 
+ 00 , we have 

d(SOCP) = u(P-SDPi) = u(D-SDP). (17) 

Moreover, the problem (UQ,) enjoys the strong duality under the further as¬ 
sumption ©, i.e., w(UQ + ) = u(D-SDP). 

Proof. As proved at beginning of Theorem[l] if u(UQ + ) < +oo, then n(SOCP) < 
+ 00 . Replacing © with the equivalent convex quadratic constraint x T Qx < t, 
we can reformulate (SOCP) as a convex QCQP, denoted by (CQCQP). Then 
we have 


u(SOCP) = v(CQCQP). 
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It is not difficult to verify that the Lagrangian dual of (CQCQP) is exactly the 
same as that of (UQ ,), i.e., (D-SDP). Then, the Slater condition guarantees 
that (CQCQP) itself enjoys the strong duality, i.e., 

w(CQCQP) = 'u(D-SDP), 

see for example, Proposition 6.5.6 in [ 8 ]. Since the Slater condition holds for 
(P-SDPi) and its conic dual (D-SDP), strong duality holds for (P-SDPi) and 
(D-SDP) [31], i.e., 

■u(P-SDPi) = u(D-SDP). 

Therefore, we obtain the equalites m- On the other hand, since linear trans¬ 
formation does not change the strong duality, it follows from Theorem Q] that 

v(UQ + ) = tj(SOCP) 

under the further assumption (JSJ . Thus, the proof of the second part is com¬ 
plete according to (El). □ 


3 Sufficient Conditions for Hidden Convexity of the General 
(QCQP) 

Now we extend the above convexity approach for the uniform quadratic opti¬ 
mization problem (UQ-|_) to the general quadratic constrained quadratic pro¬ 
gramming (QCQP): 

(QCQP) min go(x) (or ma x—g 0 (x)) (18) 

s. t. k < g.i(x ) < Ui, i=l,...,p, (19) 

where the functions gt(x ) : R™ — > R (i = 0, 1 ,... ,p) are given by 

m 

gi(x) = ^2 aijX T QjX + 2 bf x + Cj, (20) 

i =i 

with a K j e {—1,0,1}, Qj = Qj ^ 0, bi € R™ and Ci £ R, for i = 0,1, ... ,p and 
j = 1 ,..., m. (QCQP) is the simplest type of nonconvex nonlinear program¬ 
ming. It has many applications, see m and references therein. In this section, 
we establish a new sufficient condition under which (QCQP) is equivalent to 
its convex relaxation. Our new result not only generalizes a few existing re¬ 
sults for some special (QCQP) but also partially gives answers to a few open 
questions in literature. 
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3.1 One-Sided QCQP 


In this subsection, we study the one-sided (QCQP), denoted by (QCQPi), 
which is a relatively easy case of (QCQP) with h = —oo for i = 1,... ,p. Our 
first main result is as follows. 

Theorem 3 Suppose 


max 

jeJ 



span{ 6 i,..., b p } (J.A f(Qj) (J H(Qi) 




< n — 1 , 


( 21 ) 


the problem (QCQPi) is equivalent to the following convex QCQP relaxation: 


(CR) 


min 


x T QjX + ^2 a 0 jtj + 26q x + co 
j&J 


s. t. 


y; x T QjX + ^2 aijtj + 2 bjx + d < 0, 
j#J j£J 



< tj + 1 
“ 2 


j e J , 


= !,•••, P, 


( 22 ) 


in the sense that x* globally solves (QCQPi) if and only if {x* ,t*{j £ J)) := 
(x*,x* T QjX*(j £ J )) globally solves (CR), under the assumption v(CR) > 
—oo, where 


j = Ui? e i 1 ’---’™} i 


there is an i £ {0 ,... ,p} such that aij 


-!}• 


Proof. Let ( x*,t*(j £ Jj) be an optimal solution of (CR) since u(CR) > —oo 
and the feasible region of (CR) is closed. Note that K21 ) is equivalent to 

x T QjX < tj. 

Suppose x T QjX = tj for all j £ J, the proof is complete. Now, we assume that 
there is a jo € J such that x* T Qj 0 x* < t* Q . 

Let (•)' L be the orthogonal complement of the subspace (•). Notice that 
J\f ± (Q) = Tl{Q) and (A U B )-*- = A 1 - Cl B 1 . Then, (Ell) implies that 


min 

jeJ 



(span{ 6 i,..., 6 p }) J *P)^(Q i ) QA/XQ*) 




> 1 . 


Therefore, there is an xq £ (span{b 1; ..., ftp}) -1 f) 7Z(Qj 0 ) f)i^j 0 A/"(Q») and 
Xq ^ 0. Since Qj 0 >r 0, we must have x]^Qj 0 x o > 0. Otherwise, x^Qj 0 x o = 0 
and then Qj 0 x o = 0. It follows that Xq £ 1Z(Qj 0 ) f) Af(Qj 0 ) = {0}, which is 
contradiction. According to the fact x* T Qj 0 x* < t* o , the following equation 


{x* + ax 0 ) T Qj o (x* + axo) = t* 0 
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has a solution, denoted by a o- Define 

x := x* + ctoXo. 

We can see that ( x,t*{j £ J)) is feasible for (QCQPi) with the same objective 
value as ( x*,t*(j £ J)). That is, ( x,t*(j £ J)) is also an optimal solution of 
(QCQPi). The proof is complete by noting that x r Qj 0 x = t* Q . □ 

Next, we show that Theorem [3] has many applications. 

We first consider the well-known trust region subproblem 

(TRS) min x T Ax + 2 b T x 
s. t. ||.t|| 2 < 1 , 

where A = A T £ R nx " and b £ R n . If A y 0, (TRS) is already a convex 
programming problem. Otherwise, we have: 

Corollary 1 Suppose A ^ 0, the nonconvex (TRS) is equivalent to the fol¬ 
lowing ball constrained convex quadratic programming relaxation: 

min x t (A - A mi „(A)I n )a; + 2 b T x + A min (A) 
s. t. ||x || 2 < 1 . 

Proof. We can recast (TRS) as the following formulation: 

min x T (A — A m i n (A)/n)£ + 2 b T x + A m i n (^4):r T :r 
s. t. x T x < 1 . 

According to Theorem (3J it is equivalent to the following relaxation 

(TRSR) min x T (A — A m ; n (A)I n )x + 2 b T x + A m i n (A)t 
s. t. t < 1 , x T x < t , 

since u(TRSR) > —oo and the assumption KT1) reduces to 
dim {1Z{A - A min (A)/ n )} < n - 1, 

which is always true. As we have assumed A ^ 0, it holds that A m ; n (A) < 0. 
Thus, for any optimal solution of the above problem, denoted by (a;*,t*), we 
always have t* = 1 , which completes the proof. □ 

As a further extension, we study (TRS) with additional linear inequality 
constraints Hums]: 

(ETRS) min^gRn x T Ax + a T x 

s. t. ||x — Xq || 2 < u, 

bfx < Pi, i = 1 ,... ,p, 

where A = A T £ R nx ™, a , bi, xq £ R n and u, Pi £ R.. Again, we assume 
A 0 , i.e., A m i n (A) < 0. 
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Corollary 2 Suppose A >£_ 0 and 

dim |span{6 i,..., b p } 0 1Z(A — A min (A)/ n )| < n — 1, (23) 

(ETRS) is equivalent to the following linear and ball constrained convex quadratic 
programming relaxation: 

min x t (A - X min (A)I n )x + a T x + X min (A)(u + 2x T x 0 - ||a:o|| 2 ) 
s. t. \\x — cco || 2 < u, 

bfx</3i , i = l,...,p. 

Proof. Similar to the proof of Corollary[T| it is sufficient to show that (ETRS) 
is equivalent to 

(ER) min x T (A - X min (A)I n )x + a T x + A min (A)f 
s. t . t — 2 x t Xq + ||^o|| 2 < u, 
bfx< fa, i = l,... ,p, 
x T x < t. 

Firstly, it is easy to verify that the assumption (P23l) coincides with (12T1) . To 
show u(ER) > — oo, it is sufficient to prove the feasible region of (ER) is 
bounded. According to the Cauchy-Schwartz inequality, for any feasible solu¬ 
tion of (ER), denoted by (x,t), we have 

INI 2 <t<u - ||a:o|| 2 + 2x t x 0 <u - HxoII 2 + 2||x 0 ||||x||. 

Therefore, both ||x|| and t are bounded. Now, the equivalence between (ETRS) 
and (ER) directly follows from Theorem [3] □ 

Remark 1 The assumption (E3ll was first proposed by Hsia and Sheu [IS. to 
guarantee that u(ETRS) is equal to the optimal objective value of the pri¬ 
mal SDP relaxation of (ETRS). It improved the following assumption due to 
Jeyakumar and Li m 

dim {A r(A -^min (A)I n )} > dim {span {&i,..., b p }} + 1. 

Now we focus on the weighted max-min dispersion problem din]: 
max min ujA\x — zA\ , 

xeo t=i,...,p 

where 12 C R" is closed, z\,...,z p in R" are given points and u>\,... ,oo p 
are positive weights. When Q is a ball with a radius r 0 centering at Xq , the 
max-min dispersion problem can be reformulated as a special case of (QCQP): 

(WD) max s 

s. t. s < Wi|| a; - Zi || 2 , i = l,...,p, 

||a,'-x 0 || 2 < r'l 
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An approximation algorithm for (WD) was proposed in 1161117] based on the 
following SDP relaxation: 

(WD-SDP) max s 


s. t. s < U>i 






( Y x 
[x T 1 


, i = l,...,p, (24) 


I —Xq 

-xl Ikoll 2 


Y x 

V T 1 


ho, 


Y x 
r T 1 


< r Q , 


(25) 

(26) 


where Y is relaxed from xx T . It was raised as an open question in mm that 
under what conditions the SDP relaxation of (WD) is tight, i.e., u(WD) = 
v(WD-SDP). Here, we can quickly give an answer based on Theorem [2] 

Corollary 3 Suppose 


rank [zi — xq, ..., z p — .To] < n — 1, (27) 

(WD) is equivalent to the following linear and ball constrained convex quadratic 
programming relaxation: 

(WR) max s 

s. t . s < Ui(r% - 2 (zi - t 0 ) T y + ||zj - t 0 || 2 ), i = 1,... ,p, (28) 

M < ro- (29) 

Moverover, under the assumption int(f2) ^ 0, i.e., there exists anx such that 
||x — To|| < ro, we always have r(WR) = t>(WD-SDP). 

Proof. Introducing y := x — To, we see that (WD) is equivalent to 

(WD') max s 

s- t. s < Ui\\y - (Zi - Xq ) 11 2 , i = l,...,p, 

\\y\\ 2 <rl 

For (WD'), it is trivial to verify that Assumption (l2lTl reduces to (1271) . Ac¬ 
cording to Theorem [3] (WD') is equivalent to the following convex relaxation 

(WR') max s 

S. t. s < uji(t - 2 (Zi - xo) T y + ||zj - t 0 || 2 ), i = l,...,p, 
t < ro, 
y T y < t, 

by noting that the feasible region of (WR') is clearly bounded. Let ( y,t,s) 
be any feasible solution of (WR'). (j/,rg,s) remains feasible since Wi > 0 for 
i = 1,... ,p. Therefore, we always have w(WR') = v(WR). 
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The assumption int(L?) / 0 implies that the Slater condition holds for 
(WR'). Since u(WR') < +oo, according to Proposition 6.5.6 in [B], strong 
duality holds for (WR') and its Lagrangian dual. 

Consider the SDP relaxation of (WR'): 


(WD-SDP') 


max s 


S. t. S < UJi 



xo ~ Zi \ 

Wzi-xofJ 




h 0. 


< r 


2 

0 ’ 





Let ( Z,y,s ) be any feasible solution of (WD-SDP'). It follows from Tr (Z) = 
I • Z < Tq and Z — yy T y 0 that ||j/|| is bounded since 

y T y < Tr (Z) < r%. 


Then, according to the first constraint of (WD-SDP'), s is bounded from above. 
Note that Z y 0 implies that 


Za >0, i — 1,..., n, 

ZiiZjj > j — 1; • • • i fl- 

That is, each entry of Z is bounded. Then the feasible region of (WD-SDP') 
is compact. Under the assumption int(D) / 0, the Slater condition holds for 
(WD-SDP'). Therefore, strong duality holds for (WD-SDP') and its conic dual, 
see for example [51] . 

It is not difficult to verify that the Lagrangian dual of (WR') is the same 
as the conic dual of (WD-SDP'). We have w(WD-SDP') = v(WR'). 

Let (Y, x, s) is a feasible solution of (WD-SDP). We can verify that (Z, y, s ) := 
(Y — xqx t — xxq + xoXq, x — xq, s) is a feasible solution of (WD-SDP'). On 
the other hand, for any feasible solution (WD-SDP'), say, (Z, y, s), (Y, x, s) := 
(Z + xoy T + yxQ +xoXQ,y + Xo, s)) is a feasible solution of (WD-SDP). There¬ 
fore, c(WD-SDP') = v(WD-SDP). 

As a conclusion, we obtain 

u(WR) = d(WR') = t;(WD-SDP') = v(WD-SDP). 


□ 


3.2 Two-Sided QCQP 

Now we focus on the general two-sided (QCQP) ([T5[) - ([T$p . Our second main 
result is as follows. 
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Theorem 4 Suppose 


max 

jeK 


dim < span{ 6 i,... ,b p } (J U(Qi) 


< n — 1 , 




(30) 


the problem (QCQP) is equivalent to the following linear and second-order 
cone constrained convex quadratic programming relaxation: 


(CR') min x T QjX + aojtj + 2b^x + Cq 

j?K j£K 

s. t. k < Qijtj + 2bfx + Ci < Ui, i = 1 ,... ,p, 
jeK 

tj + l . 


Qj x 

tp-1 


< 


j e J, j g K, 


in the sense that x* globally solves (QCQP) if and only if ( x*,t*(j € K)) := 
(x*, x* T QjX*{j £ K )) globally solves (CR'), under the assumption v(CR') > 
—oo, where 


K = , m} | aoj = —1 or ajj / 0 for some i € { 1 ,... , p}} . 

Our first application of Theorem [4] is for the general problem (UQ) rather 
than (UQ+). We first consider the positive semidefinite case. 

Corollary 4 Suppose Q >: 0 and 


rank [bi,... ,b p ] < rank(Q) — 1 , 


The problem (UQ) is equivalent to (SOCP) in the sense that x* globally 
solves (UQ) if and only if(x*,t*) := (x*,x* T x*) globally solves (SOCP) (0). 

Then we study the indefinite case of (UQ). Let Q = Q± — Q 2 be a decomposition 
such that Qi >z 0, Q 2 h 0, rank(Qi) = ri, and rank(Q 2 ) = r 2 . We can 
reformulate the indefinite (UQ) as: 

max x T Q\x — x T Q 2 X + 2b^x + do 

s. t. k < x t Q\x — x T Q 2 X + 2 bjx < Ui, i = 1 ,... ,p. 

Applying Theorem 2] yields the following result. 

Corollary 5 Suppose 


rank [ 61 ,..., b p \ < min{n, r 2 } — 1 
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the indefinite (UQ) is equivalent to the following second-order cone program¬ 
ming relaxation: 


(SOCP') max l£ t« it — £2 + a; + d 0 

s. t. li < ti — t 2 + 2 bfx < Ui, i = 1 ,... ,p, 
t\ + 1 



in the sense that x* globally solves (UQ) if and only if 

(x*,tl,t* 2 ) := (x*,x* T Q 1 x*,x* T Q 2 x*) 


globally solves (SOCP 7 ). 

The two-sided ball constrained trust region subproblem was first proposed by 
Stern and Wolkowicz EH], see also U: 

(TTRS) min f(x) = ^ x T Ax + b T x 
s.t. a < x T x < /3, 

where A = A T G M nxn , b £ R n and —00 < a < /3 < + 00 . Directly applying 
Theorem 0] we have: 

Corollary 6 The problem (TTRS) is equivalent to the following linear and 
second-order cone constrained convex quadratic programming relaxation: 

(CTTRS) min -x T (A — A m i n (A)/)x + b T x + X m i n (A)t 
s.t. a < t < (3, 



in the sense thatx* globally solves (TTRS) if and only if (x* ,t*) := (x*, x* T x*) 
globally solves (CTTRS). 

Very recently, Bienstock and Michalka [3] studied the following variants of 
trust region subproblem: 

(VTRS) min x T Qx + c T x 

s. t. ||x < n, i G I, 

Ik-Mjll > rj, j G J, 

X G P, 

where P = {x | a^x <bk, k = 1,..., m}. It was showed in |3] that (VTRS) is 
polynomially solvable as long as the number of faces of P within the ellipsoids 
is polynomial. 
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Decomposing Q = (Q - A m i n (Q)I n ) + A m ; n {Q)I n and then applying The¬ 
orem U yields the following new sufficient condition to guarantee the hidden 
convexity of (VTRS): 

Corollary 7 Suppose 

dim jspan{ai,.. .,a m ,pi(i G G J)} (JA1(Q - A min (Q)/ n )| < n - 1, 

the problem (VTRS) is equivalent to the following linear and second-order cone 
constrained convex quadratic programming relaxation: 

min X (() A m i n(Q)In'jx T C X -t- A min(Q)t 

S. t. t - 2pjx + \\ni\\ 2 <Ti, i G I, 
t~ 2pjx + \\ixj\\ 2 > rj, j G J, 



a%x < b k , k = 1 ,..., to , 


to t/ie sense thatx* globally solves (VTRS) if and only if(x*,t*) := (x*, x* T x*) 
globally solves the convex relaxation problem. 

Finally, we consider the extended p-regularized subproblem (p-RS m ) with 
to. inequality constraints HBj: 

(p-RS m ) min xeH n \x T Qx + c T x + -||a;|| p 

2 p 

s. t. li < afx + bi < Ui, i = 


where Q = Q T G R” xn , c, G M”, bi, k, m G M, a > 0 and p > 2. Hsia et 
al. (T8j showed that if to is fixed then (p-RS m ) is polynomially solvable when 
p = 4. The complexity for the general p > 2 remains unknown, see m- 

Again, based on the matrix decomposition Q = (Q—X m i n (Q)I n )+X m i n (Q)I n , 
we apply Theorem[4]to establish a new polynomially solvable class of (p-RS m ) 
for p > 2: 


Corollary 8 Under the condition that 

dim |span{ai,..., a m ) R.{C) A m j n j- n 1, 

the problem fp-RS m j is equivalent to the following linear and second-order 
cone constrained convex quadratic programming relaxation: 

min -x t (Q - A mi n (Q)I n )x + c T x + -t% + A min t 

2 p 

s. t. li < af x + bi < Ui, i = 1,... , to, 


x 

t -1 
2 


< 


t "j- 1 


in the sense thatx* globally solves (p-RS m ) if and only if (x* ,t*) := (x* ,x* T x*) 
globally solves the convex relaxation problem. 
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4 Approximation Algorithm for Convex-Constrained (UQ_|_) 

As shown in Theorem [2] u(UQ + ) = u(SOCP) under Assumption ©. Suppose 
© does not hold, it is not difficult to find examples satisfying u(UQ + ) < 
v(SOCP), see [3]. Thus, it is natural to ask what is the quality of v(SOCP). In 
this section, we answer this question for the convex-constrained case of (UQ+) 
where li = — oo, i = 1,... ,p. More precisely, a new approximation algorithm 
based on (SOCP) is proposed for solving (UQ+). It should be noted that the 
approximation algorithms for the general nonconvex quadratic optimization 
problem with ellipsoid constraints n work well for (UQ+). However, as 
we see below, for (UQ+), the existing approximation bounds can be greatly 
improved. Moreover, the existing approximation algorithms j3T)ll20l are based 
on SDP relaxation, while our new algorithm is based on SOCP relaxation. 

Theorem 5 Under the assumption that the origin 0 is in the interior of the 
feasible region of (UQ+), do = 0 , U = — oo and Ui < +oo for i = 1,... ,p, we 
can find a feasible solution x in polynomial time satisfying 

u(SOCP) > u(UQ+) > f 0 {x) > ( l ~ 1 ^ • u(SOCP) > 0, (31) 

V v 2 + 7 / 

where 7 = max i= i p , ^ 

'Jui-di + WQ-TibiP 

Proof. The assumptions 0 is in the interior of the feasible region and do = 0 
imply that u(SOCP) > u(UQ+) > 0. Under the assumption m < +00 for 
i = 1,... ,p, the feasible region of (SOCP) is compact and hence u(SOCP) < 
+ 00 . Let ( x*,t*) be an optimal solution of (SOCP). It is sufficient to assume 
x* t Qx* < t* since otherwise, v(UQ + ) = f(SOCP). Now we have 

t* + 2 b^x* = u(SOCP), t* + 2bfx* < Ui - di, x* t Qx* <t*, i = l,...,p. 

Let y g M" satisfy 

x* t Qx* + y T Qy = t*. 

Then we have y 0, y T Qy > 0 and 

(a 2 + 1) (x* t Qx* + y T Qy + 2b^x*) = (a 2 + l) • w(SOCP). (32) 

Since x* T Qx* + 2b^x* < n(SOCP) and y T Qy > 0, there exits a real value 
a > 0 such that 

(x* + ay) T Q(x* + ay) + 26q (x* + ay) = u(SOCP). (33) 

The equation (15^1) minus the equation (1551) equals 

(ax* — y) T Q(ax* — y) + 2ab^(ax* — y) = a 2 ■ r;(SOCP). (34) 

x* + ay ax* — y 

\/l + a 2 ’ y/1 + a 2 


Define 


si = 
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Then (TTTT1) and (ETH) can be recast as 

2b% Sl _ n(SOCP) 


si Qs i + 


s 2 Qs 2 + 


2 1 + a 2 ’ 

T „ „,2„ 


\/l + a 


2ab£s 2 _ a 2 v(SOCP) 


\/l + a 2 


Since 


sf Qs\ + s 2 Qs 2 = t*, x* = 


1 + a 2 

si as 2 


(35) 

(36) 


\/\ + a 2 \/l + a 2 ’ 


it follows from t* + 2 bj x* < Ui — d,; that 


26/ Si T< „ 2abl s 2 , 

Si Qsi H , + s 2 (js 2 ^ 0 < it* — dj, 4 — 1,... ,p, 


Vl + a 2 


or equivalently, 




\/l + a 2 


tt i 0:0 ^ 

( ? 2S2 + lfe^ 


Wi-di + ||Q 2 6i|| 2 Ui-dj + IIQ ^6i|| 2 
Then we have 


< 1, i = l,...,p. 


(1 + a 2 ) 


mm < max 


Q~*bi 








max 


qQ 2 bj 


Q-S 2 + ^ 


*=1’->p Wj-dj + IIQ 2 6j|| 2 i=i,...,p Ui-di + ||Q 2 ^ll 2 

. f n lT a 2 
< mm < 1 + a , --— 


1 + a 2 , if |a| < 1, 
1 + -K , otherwise, 


< 2 . 


Therefore, there is an index j £ {1,2} such that 

WQisj/tj + Q-ibiW 


where t2 = Tife? 


- di + ||Q 2&,j||2 
. Define 


< "n/ 2 , *= l,...,p. 


(37) 


£ := 


— J O' 3 


s-/t -, if&o s yAj-° 


— Sj/tj, otherwise, 
r := max re[0 ,lj: fi(jx) < u i: i = l,...,p 


I r rn n \\tQzx + Q 26 i|| 1 . t 

= max < r € [0,1J : — -j = = - = <1, 4 = 1,... ,p 


\Jui-di+\\Q 2 bi 11 2 
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According to (IH71) , it holds that 


WQzx + Q ibj |[ 

\Jui-di + \\Q~hi\\ 2 

+ Q~^ b i\\ II - {Q^sj/tj + Q-ibi) + 2Q~hi\\ 


< max ■_ 

Ui - (k + WQ-^biW 2 \jui- di+ \\Q-^bi \\ 2 

< max L/9 V2+ 2| 1Q~^H - 

[ yu* - di + 

, y /o , 2 ||Q-*M| 

lui- di + \\Q~^bi \\ 2 


Therefore, for any r G [0,1], we obtain 

\\tQ^x + Q~ibj\\ 

\Jui- di + \\Q~hbiW 2 

\\t{Q^x + Q-^bi) + (1 - t)Q~?1 


yJui-di + \\Q 2&i|| : 


< t [ y/2 + 


nQ-*bi\\ 


Ui-di + \\Q *bi || 2 , 


+ (l-r) 


ll<T*M 


b-rfi + ||Q 2 11 2 


Since the origin 0 is in the interior of the feasible region of (UQ+), we have 


Ui — di > 0 for i = 1,... ,p. It follows that 


i = 1 ,... ,p, we have 


\j Ui — 


11QII&ill 


< 1. Now, for 


di+IIQ-sMI 2 


HtQU + Q 26 ^ 


< 1 , 


/ui-di + ||Q 2&i|| 2 
or equivalently, fi(fx ) < Ui, as long as 


r < 1 


|| Q~*bi 



V2 + 


IIQ-^ill 


\Jui-di + \\Q 2&i|| 2 / ' \ ^/ui-dj + HQ zbiW 2 
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According to the definition of r, we obtain 


T > 



Ui-di + \\Q 2 &i || 2 



where the equality holds as /i(y) = (1 — 7)/(-\/2 + 7 ) is a decreasing function 
for 7 G [0,1). 


Now, we conclude that 


fo(rx) =t 2 x t Qx + 2rb^x 


>t 2 x t Qx + 2 r 2 b^x 
>t 2 x t Qx + 2t 2 6q sj/tj 
=T 2 {s-fQs-j+ 2tjb%sj)/tj 2 
=t 2 ■ v(SOCP), 


(38) 

(39) 


(40) 


where the inequality (1551) follows since b^x > 0 and t > r 2 (r G [0,1]), the 
inequality (1551) holds as b^x > Sj/tj, and the equality (15(1 follows from (1551) 
and (l36l) . □ 

Remark 2 To satisfy the assumption that 0 is in the interior of the feasible 
region of (UQ + ), it is sufficient to find any one of the interior points of (UQ+) 
and then translate the origin there. 


5 Approximation Algorithm for Finding the Chebyshev Center of 
the Intersection of Balls 

As an application of (UQ+), Beck [51[5] studied the problem finding the Cheby¬ 
shev center of the intersection of given balls, i.e., the min-max problem (CC) 
O. Replacing the inner maximization problem with its Lagrangian dual, which 
is the following minimization SDP problem: 


(SDP(z)) min t 



A * G M+, i = l,...,p, 
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Beck [3] reduced the min-max problem (CC) to a double minimization prob¬ 
lem, which turns out to be convex programming problem: 


(DCC) 


min d(SDP(z)) + ||z || 2 

Z 

= min t + ||z || 2 

t,\,Z 

/ (—1 + EiLl ~ — Ef=l ^i a i 

l t + ELiMIM 2 - 

A i e R+, i = 1 ,... ,p. 



ho, 


Moreover, by noting that the optimal objective value of above convex SDP 
problem can be attained at a solution (t*,X*,z*) satisfying Eo=i A* = 1 and 
z* = Ei=i A* a,;. Beck j3] showed that (DCC) is further equivalent to the 
following standard convex quadratic programming problem: 


min A J2 Xi ( r i “ INI 2 ) 

i =1 
V 



2 


s. t. y> = l, \ > 0 , i = 1 ,... ,p. 

i —1 


Based on the strong duality theory for (UQ+), Beck [4] showed w(CC) = 
f;(DCC) when p < n. For the hard case p > n, solving (DCC) always yields a 
heuristic Chebyshev center z. To our knowledge, it remains unknown what is 
the quality of r>(DCC) and z. In this section, we establish the first approxima¬ 
tion ratio between v(CC) and u(DCC). Moreover, the quality of the solution 
z returned by Beck’s approach is also guaranteed. 

Theorem 6 Under the assumption that int(C) ^ 0, we can find a solution z 
in polynomial time satisfying 


i>(DCC) > max ||ar 


z || 2 > u(CC) > 



2 

■v(DCC), 


(41) 


where 7 is equal to the optimal objective value of the following convex program¬ 
ming problem: 


min max 

xGM 71 i=l,...,p 


II*-«»ll 


Ti 


(42) 


Moreover, let d max = maxj^i,,..^ ||a* — <Zj|| and r m ; n 


min^i,...^ rj, we have 


7 < 


2 (n + 1 ) r, 


dmax dmax 


V2 


(43) 
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Proof. According to Theorem EJ u(SDP(A)) = u(SOCP(^)), where SOCP(^) 
is the second-order cone programming relaxation proposed in Section 2: 


SOCPO) 


max t — 2z t x 

s. t . t — 2 afx + |||| 2 < r 2 , 



= P, 


Suppose 0 £ int(l7), according to Theorem [5j for any z £ I", we have 
u(SDP(z)) > max{||a;|| 2 — 2 x T z} > r 2 • x(SDP(z)), 

x£f2 


where 

1 - max, = i r i ! p ^ 

\/2 + maxj=i. . ^ 

Since 0 £ int(l7), we have ||aj|| < r* for * = 1,. .. ,p. That is, r > 0. Notice 
that it is trivial to see r < 1. Therefore, t 2 < 1 and then it holds that 

min{?;(SDP(z)) + ||z|| 2 } > minmax||x — z\\ 2 > r 2 min{w(SDP(z)) + ||z|| 2 }, 

2 : z x^lQ z 

or equivalently, 

z;(DCC) > u(CC) > r 2 • u(DCC). (44) 

According to the definition of (SOCP(z)), we have 

v(SOCP(z)) > ||x || 2 - 2 z T x, V x £ n. 

Let z be an optimal solution of (DCC). Since x(SDP(z)) = t;(SOCP(z)), it 
holds that 

u(SDP(z)) + ||z || 2 > ||x — z|| 2 , V x £ 17. (45) 

Therefore, combining (TH1) and fi5l) yields 

i;(DCC) > max ||s - z || 2 > v(CC) > r 2 • v(DCC). (46) 

x£f2 

Let xo be any interior point of 17, i.e., ||xo — a*|| < r*, i = 1 ,.. . ,p. Define 
(CC') min max \\x — (z — Xo)\\ 2 , 

z x£f2(x o ) 

where J7(xo) = {x £ R" | ||x — (a* — xo )|| 2 < r 2 , i = 1,... ,p}. Then it is 
trivial to see that 0 £ int(f7(xo)) and u(CC) = u(CC'). 

We write the second-order cone programming relaxation for the inner max¬ 
imization problem of (CC 7 ) as follows: 

SOCP'(z) max t — 2{z — Xq) t x 

s. t. t- 2 (oj - x 0 ) T x + || aj - x 0 || 2 < r 2 , i = l,...,p, 
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Let (x*,t*) be an optimal solution of S0CP( 2 ). Define 

x := x* — xo, t:=t* — 2x^x* + x^Xq. 

We can verify that (x, t) is a feasible solution of S0CP'( 2 ) and hence 

v(SOCP'(z)) > t- 2(z - x 0 ) T x = v(SOCP(z)) + ||z || 2 - \\z - x 0 \\ 2 . (47) 

On the other hand, let (x*,t*) be an optimal solution of S0CP'( 2 ). Define 
x := x* + xo, t := t* + 2xqX* — XqXq- 
We can verify that (x,t) is a feasible solution of SOCP(z) and hence 

v(SOCP(*)) > t - 2 z T x = u(SOCP '(*)) - || 2 || 2 + \\z - x 0 f. (48) 
Combining (TX71) and (H51) yields 

w(S0CP( 2 )) + || 2 || 2 = v(S0CP'( 2 )) + \\z ~ xof- (49) 

Denote by SDP'( 2 ) the SDP relaxation for the inner maximization problem of 
(CC'). According to Theorem [ 2 ] we have u(SDP , (z)) = i>(S0CP'( 2 )). Then, 
the equation (PI becomes 

u(SDP( 2 )) + || 2 || 2 = t,(SDP'(*)) + \\z - soil 2 - (50) 

Similar to (DCC), we define 

(DCC') min u(SDP'( 2 )) + || 2 -a: 0 || 2 . 

Z 

Then the equation (15U1) implies that 

n(DCC') = v(DCC). 

Since 0 £ int(f2(xo)), according to flU . we have 

i>(DCC') > _max \\x — (z — £o )|| 2 > 'u(CC') > t(x o ) 2 • i>(DCC'), 

x(zf2(x o) 

or equivalently, 

u(DCC) > max \\x — 2 || 2 > v(CC) > t(x o ) 2 • u(DCC), (51) 

x£f2 

where z is an optimal solution of (DCC') and 

1 — maxj = i „ 

Since the inequality m holds for any xo £ int (.!?), we can choose xq to 
maximize the lower bound t(xq) in int ( 12 ), that is, 



t(x o ) 2 = 


2 
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where 


7 = 


inf max 

xo£int(f2) i=l,...,p 


ll-'CQ - ai|| 
n 


= min max 

xo£f2 


Iko - aj\\ 
n 


(52) 


Since ||xo — a* || < r,; for any x G fi, we have 7 < 1. Let x* be an optimal 
solution of ( 021 ) • Notice that the optimal objective value of (B21) provides a 
lower bound of the problem (|521) . Then, we have maxj = i v .. iP ^ ~ ai ^ < 1, that 
is, x* £ fl. Therefore, the problems flUD and (El are actually equivalent. 

According to the definition of r m i n , we obtain 


7 < min max 


11^0 - Qi|| 

Train 


-- min max ||xo — a,i\\. 

T min 


Since the objective function of the above inner maximization problem is con¬ 
cave with respect to a*, we have 


7 < -- min max llx’o 

T min ®oeR n aGC 


a|| 


where C is the convex hull of the points a± ,...,a p . According to Example 
3.3.6 in [IT], we have 


min max ||xo 

i 0 ei n aec 


< d. 


— u max 


2(n+l) 


and then the proof is complete. 


□ 


6 Conclusion 

In this paper, we have thoroughly studied the uniform quadratic optimization 
problem (UQ+), which is a special quadratic constrained quadratic program¬ 
ming (QCQP) sharing the same positive definite Hessian matrix Q. We show 
that (UQ+) is NP-hard when p > n + 2, where n and p are the numbers 
of variables and constraints, respectively. We proposed a second-order cone 
programming (SOCP) relaxation for (UQ+) and showed that it is tight un¬ 
der a new assumption, which is slightly weaker than the existing assumption 
p < n. Now the complexity of (UQ+) remains unknown only when p = n + 1. 
Then, we extended the SOCP relaxation approach to (QCQP) and established 
a new sufficient condition under which (QCQP) is hidden convex. As further 
applications, this new condition not only generalized some well-known results 
(for example, the trust region subproblem and extensions, (UQ) with positive 
semidefinite Q and indefinite Q), but also partially gave answers to a few open 
questions (for example, the max-min dispersion problem and the extended p- 
regularized subproblem). For convex constrained (UQ+), we proposed an im- 
roved approximation algorithm based on our SOCP relaxation. Moreover, the 
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new approximation bound is independent of p and n. As a further application, 
we succeeded in establishing the first approximation analysis for Beck’s convex 
approach to find the Chebyshev center of the intersection of p balls. Again, 
the approximation bound dependents only on the distribution of the given 
balls rather than p and n. The future work may include more applications of 
our sufficient condition for general (QCQP) and further improvement of our 
approximation analysis. 
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